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Abstract 

The nucleon-nucleon (NN) t-matrix is calculated directly as function of two 
vector momenta for different realistic NN potentials. To facilitate this a for- 
malism is developed for solving the two-nucleon Lippmann-Schwinger equa- 
tion in momentum space without employing a partial wave decomposition. 
The total spin is treated in a helicity representation. Two different realistic 
NN interactions, one defined in momentum space and one in coordinate space, 
are presented in a form suited for this formulation. The angular and momen- 
tum dependence of the full amplitude is studied and displayed. A partial wave 
decomposition of the full amplitude it carried out to compare the presented 
results with the well known phase shifts provided by those interactions. 



PACS number(s): 21.45.+V, 13.75.Cs 



1 



Typeset using REVdp^X 



I. INTRODUCTION 



Experience in three- and four-nucleon calculations [|I|-|3| shows that the standard treat- 
ment based on a partial wave projected momentum space basis is quite successful but also 
rather tedious, since each building block requires extended algebra. For example, the rep- 
resentation of the various permutation operators (i.e. the transformations among angular 
momenta belonging to different Jacobi momenta) is very involved |4[] and requires intri- 
cate numerical realizations ||. On the other hand, we demonstrated in the relative 
ease with which a 3-boson bound state could be calculated in the Faddeev scheme avoiding 
an angular momentum decomposition altogether. Instead of solving a large set of coupled 
two-dimensional integral equations in the standard partial wave framework, in a three di- 
mensional (3D) formulation only one single integral equation in three dimensions had to be 
solved. We had the same positive experience with the scattering of three bosons 0, where 
the algebraic formulation and numerical implementation of the one Faddeev-equation was 
much simpler than the machinery of the approach based on partial wave decomposition. 

The input to these three-body calculations without angular momentum decomposition 
are two-body t-matrices, which are off the energy shell. In addition to the off-shell energy 
they depend on the magnitudes of the initial and final momenta and the angle between the 
two momenta. In || we showed that the two-body t-matrices can be obtained very eas- 
ily by solving a two-dimensional Lippmann-Schwinger (L.S.) equation instead of preparing 
and handling quite a few one- dimensional L.S. equations for each angular momentum state 
separately. The off-shell energies required in a three-body calculation lie between the total 
three-body energy and minus infinity. Interestingly we found in || that the t-matrix at large 
negative energies is nearly equal to the real part of the t-matrix at corresponding positive 
energy. This insight, which was also pointed out in ||, explains directly that for bound 
state calculations a surprisingly large number of partial waves is necessary for convergence. 
In addition, recent three-nucleon (3N) Faddeev calculations |10| based on realistic nucleon- 
nucleon (NN) forces and performed between 100-200 MeV nucleon laboratory energy reach 
the very limits of present day computational resources, since the number of angular momen- 
tum states involved increases dramatically. Regarding all the above considerations it is a 
challenge to incorporate spin into the previous calculations |6|-§[ and perform few-nucleon 
calculations without angular momentum decomposition. As an aside, in the Monte Carlo 
calculations in configuration space ]TT[ this is anyhow a standard procedure. 

In the case of NN scattering calculations without partial wave decomposition have already 
been realized several times [FJ,I3]. In those studies the spin states were treated in their 



individual m-representations. A different approach is based on a helicity representation [JT4 
where the helicity related to the total two-nucleon spin was introduced. In the present work 
we follow this basic philosophy, however end up with somewhat different final equations 
to be solved. We apply our formulation to the realistic NN potentials Bonn-B [103], as 



well as the Argonne AV18 JT6[ potential and check the accuracy of our calculations by 
comparing to phase shift parameters obtained from corresponding partial wave calculations. 
The potentials we chose for implementation are given in an operator form, and thus directly 
applicable to our formulation without partial wave decomposition. Other modern NN forces, 



like Nijm I and II |T7| or CD-Bonn [|T^], are parameterized for each angular momentum state 
separately and thus are not as useful for an approach without partial wave decomposition. 
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This article is structured in the following way: In Section II we present the formalism. 
The helicity representation of the potentials and the final form of the Lippmann Schwinger 
equations are displayed in Section III. The implementation of two modern realistic NN forces 
into our formulation is given in Section IV. The connection to the standard partial wave 
representation is given in Section V. Our calculations and results are presented in Section 
VI, and we conclude in Section VII. 



II. FORMULATION FOR TWO NUCLEON SCATTERING BASED ON 

HELICITIES 

First, we introduce a helicity basis for the total spin S of two nucleons. For the z-axis 
being quantization axis, the total spin state has the well known form 



\zsA) = ]T C(-^S;m 1 m 2 A) 



.1 

2— mi 
2 



z l -m 2 ) (2.1) 



Applying the rotation operator 

R(q) = e -iS^ e -iS u e^ (2.2) 

where S z , S y are components of the spin operator S = + c 2 ), leads to the general state 

\qSA) = R{q) \zSA) (2.3) 
This is eigenstate to the helicity operator S • q: 

S ■ q \qSA) = A \qSA) (2.4) 

This follows simply from R(q)S ■ 5i? _1 (g) = S • q. Next we define momentum-helicity states 
as 

|q;gSA> = |q> |gSA> , (2.5) 

where q is the relative momentum of the two nucleons. 
The parity operator P acts as 

P\q-qSA) = \-q;qSA) (2.6) 
on the momentum helicity eigenstates. Consequently parity eigenstates are given as 

|q; qSA) n = + Vn P) |q; qSA) , (2.7) 

where rj n = ±1. 

Combining Eq. (|2.7|) with two-body isospin states \tmt),we introduce antisymmetrized two- 
nucleon states as 

|q; qSA; t) m = -^(1 - A 2 )^(l + r^P) |q; qSA) \t) 

= ^(l-vA-) S+t )\t)h,qSA) w (2.8) 



To arrive at Eq. (|2.8|) we use well known properties of two-nucleon spin- and isospin states. 
In order to evaluate the normalization of the states in Eq. ( p.8|) we need the relation between 
\qSA) and |— qSA). We use the definition of Eq. ( |2.3| ) for |— qSA) and the Wigner D- function, 
namely 

D s A , A (q) = (zSA'\R(q)\zSA) 

= e-^d s A , A (6), (2.9) 

to obtain 

\-qSA) =Y,D s A , A (-q)\zSA') 

A' 

= J2 e- i{nw)A 'd s A , A {7r - 9) \zSA') 

A' 

= J2e- l ^ A '(-) s+A 'd s A ^_ A (6) \zSA') 

A' 

= (-) s E4,-a(P)I^a') 

A' 

= (-f\qS-A) (2.10) 
Now the normalization of the states given in Eq. fl2.8|) can be worked out as 



^(q'; q'S'A'; t'\q; qSA; t) m = \{l - ^,{-) s ' +t '){\ - Vn (-) S+t )S tH ^(q'; q'S'A'\^ qSA) n 

= (1 - vA-) S+t )$t'tSv n ,vJs's 

- q)5 A 'A + ^(-)^(q' + cj)5a',-a) (2.11) 



Using this result it can also be verified that the completeness relation of the states defined 
in Eq. ( |2.8j ) takes the form 

Y, I d 3 q\q;qSA;t)™ - TO <q; qSA; t\ = 1. (2.12) 



SAnt 



Equipped with the above given basis states one can formulate the Lippmann Schwinger 
integral equation. We define as matrix elements 

TfKq', q) = (q' ; q'SA'; t\ T |q; qSA; t)* a (2.13) 
VZa (q', q) = ™ (q'l q'SA>; t\ V |q; qSA; t)™ (2.14) 

Then using Eqs. ( p . 1 2j - p7T4|) the operator equation T = V + VGqT, which has a driving term 
the nucleon-nucleon (NN) potential V takes the form 

Tffiq', q) = VfM, q) + - £ / d?q" ^,(q', <f)G (<f)T&ktf , q) (2-15) 

4 A" J 

We now distinguish between the two cases for total spin ,5 = and S = 1. For S = 
Eq. (|2.15|) is one equation, similar to the one discussed in Ref. [§] for the bosonic case. For 
S = 1 there are 3 coupled equations to each A = —1,0,1. We use the property that V 
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conserves spin and isospin, which is valid to a high degree of accuracy. The coupled sets 
of equations in Eq. ( p. 15 ) can be further reduced as shown below. Using Eq. fl2.8|) and the 
parity invariance of V one obtains 

VfM, q) = ~(1 - vA-) S+t f (t\ (q; ^A'| V |q; qSA) n \t) 

= v^l - *r(-) s+ *) <t| (g'5-A'l (q'| V |q) ff |$SA) |t) , (2.16) 

where 

|q)^-^(|q>+^|-q». (2.17) 
This expression can be connected to V^f A (q', q) using Eqs. ( |2.10|) and (|2 . 1 7|) with the result 

vS? A (q',q) = vA-fvzl* (-q',q). (2.18) 

A corresponding relation is also valid for the matrix element TJ[^(q', q). Similarly one finds 

^ S - A (q' q) = vA-) s v^^, - q ). (2.19) 

We can now simplify the set of coupled equations ( [2.15|) for S = 1 in the following way: 

/ d 3 q" ^(q', q")G (/)T^ A (q", q) 

= / H'rhn 3 ^^,-^^')^-) 3 ^ (-q",q) 
= / rf 3 g"^(q / ,-q")G (g")r i ':f (-q",q) 

= / rfy y A f/(q', q") (q", q)- (2-20) 

This leads to 

T A ^(q', q) = F A ^(q', q) + i / dY t^f (q', q'OW (q", q) 

+ ^ / A" YJ? (q', q")G (g")^(q", q). (2.21) 

We now have for the case S 1 = 1 two coupled equations, namely for A' = 1,0 for each A. 
Moreover due to Eq ( |2.19p , which is also valid for the t-matrix T, it is sufficient to consider 
A = and 1. As an aside it should be mentioned that because of the relation 7r(— ) s+t = — 1 
the isospin quantum number t is fixed, once n and S are chosen. 

Finally we need to calculate the physical T-matrix element expressed in terms of the 
states 

|z/iz/ 2 m 1 m 2 q) a = -^=(1 - P 12 ) \v x v 2 m x m 2 ^i , (2.22) 

where V\,v 2 and m^mj are the magnetic isospin and spin quantum numbers, q describes 
the relative momentum of the two nucleons, and P\ 2 is the permutation operator between 
the particles 1 and 2. The physical T-matrix element is then given as 
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a (viV2m' 1 m 2 q!\T\viV2m 1 m 2 q)a = {vxP 2 m' 1 m' 2 ((\T{l - P\ 2 )\v\V 2 m\m 2 q) . (2.23) 

The above equation now has to be expressed in terms of T^*(q', q). In doing so we encounter 
the matrix element 

(q'SA' \ mi m 2 ) = Y.C(\ l -S' ] m 1 m 2 k Q ) (q'SA' \zS'A ) 
s' 11 

= C{\\s ] m 1 m 2 k )e^Ul Ao {-6 1 ) 

= C(~S; mi m 2 A )e^'dl A! (e'). (2.24) 

After some straightforward algebra one finds 
"■"(q'; q'SA'; t\piv 2 mim 2 q) a 

= - vA-) s+t )^(tWi^}A^, q'SA'\(i - P 12 )\ mi m 2 q) 

= w)C(^S; m im2 A o )e^d s Ao U0') 

(1 - V7T (-) S+t ) (%' - q) + VM + q)) • (2.25) 

We then insert the completeness relation, Eq. fl2.12|) , twice into Eq. (|2.23|) , and using 
Eq. ( p.25| ) obtain after some algebra 

a (viV2m' 1 m' 2 q'\T\i> 1 p 2 m 1 m 2 q) a 

= \e-W-^ £(1 - ^{-) s ^)C{\\t- v x v 2 f 

Snt 

C ^2~2 S] m/ i m 2 A o) C, (^ 5 '' m i m 2 A o) 

E ^(«')4A(^!(q' q)- (2.26) 

A'A 

For our calculation we choose q = z and use d AoA (0) = S AqA . Furthermore, as will be shown 
in Section III, the azimuthal dependence of q' and q can be factored out and one finds 

T$L (q', q) = e^-tiTtf* (q 1 , q, 9') (2.27) 

With this Eq. (|2.26|) can be written as (note q' = q) 

^ (v x v 2 m! x m! 2 qc[\ T \v\V 2 vn x m 2 q) a 

= -e-^~^' £ C{\\t- v x v 2 )\\ - r,^-) S+t ) 

4 Snt A A 

C(~~S; m' im ' 2 A' )C(~~S-, mi m 2 A ) £ d s A , A ,(9')T^ o (q, q, 9') (2.28) 
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III. GENERAL STRUCTURE OF THE POTENTIAL OPERATOR AND FINAL 
FORM OF THE SCATTERING EQUATION 



As it is well known, rotational-, parity-, and time reversal invariance restricts any NN 
potential V to be formed out of six terms [20f|, which are given by 



W 1 


= 1 






w 2 


= ("1 


+ tr 2 ) 


■N 






N(T 2 


■ N 




= cri ■ 


P(T 2 


P 




= <Ti • 


K<T 2 


■ k 


w 6 




P(T 2 


K + Kcr 2 - P 



(3.1) 

where the unit vectors N, P and K are defined in terms of the relative momenta q and q' 

as 

N = ^4 (3.2) 
|q x q'| 

P = (3.3) 

|q+q| 

k = ^- q | . (3.4) 

|q - q| 

The most general potential is then given as linear combination of those operators, namely 

(q'|y|q) = y(q',q) = ^(q 1 ,q,j)W i} (3.5) 

i=i 

where Vi(q', q, 7) are scalar functions depending on the magnitudes of the vectors q', q, and 
the angle between the two, 7 = q' ■ q. In Section IV we will provide specific realizations of 
the potential V(q', q). 

The evaluation of the matrix elements of Eq. (|2.14 ) for the set of operators Wi is in 



principle possible. However, it is simpler to define a different set of six independent operators, 
Qi, which is more adapt to our choice of basis states. We introduce the following set of 
operators 





= 1 






n 2 


= s 2 








= s • 


g'S 


q' 




= s- 


q'S 


Q 


n 5 


= (S 


■q? 


(S-qf 


n 6 


= s • 


qS- 


q. 



(3.6) 

The two sets of operators given in Eqs. fl3.1p and ( ^.6[ ) can be easily related via 
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ir, X-i-A. (3-7) 

j 

The transformation matrix and its inverse are explicitly given in Appendix A. 

Using the operators fij, the spin-dependent part of the matrix elements of Eq. (|2.14|) can 
then be reduced to the evaluation of matrix elements of the type (q'SA'\ f2j \qSA), namely 



(q'SA'\ 




\qSA) 


= (q'SA' \qSA) 


(q'SA'\ 


n 2 


\qSA) 


= S{S + 1) (q'SA' \qSA) 


(q'SA'\ 


n 3 


\qSA) 


= A' 2 (q'SA' \qSA) 


(q'SA'\ 


n 4 


\qSA) 


= A' A (q'SA' \qSA) 


(q'SA'\ 


fi 5 


\qSA) 


= A' 2 A 2 (q'SA' \qSA) 


(q'SA'\ 


fi 6 


\qSA) 


= A 2 (q'SA' \qSA) . 



Now we are left with determining the overlap of the states defined in Eq. ( |2.3|) . Using 
Eq. ( ^9|) this gives 

(q'SA' \qSA) = ]T (q'SA' \zSM) (zSM \qSA) 
M 

= E^f/A'(^ / 0)* J Df/ A (^0) 

M 

= £ e iM ^d s MA ,(9')d s MA (9) (3.9) 

M=-S 

A special case is the situation where q points in ^-direction, where one obtains the simple 
form 

(q'SA' \zSA) = e^'-^di^O') (3.10) 
With these preliminaries the potential matrix element of Eq. (|2.14j ) can be written as 
^(q',q) = V / 2(l-^(-) s+t ) 

E 4* (<i' SA '\ «i I^a) (*l (q'l v * Iq), l*> ( 3 - n ) 

The coefficients Aij and the matrix elements (q'| Vi\q) n depend on the angle between q 
and q' and thus on 7. Therefore their azimuthal dependence is determined by cos(0' — 
0). Furthermore, the matrix elements Qj depend on the angles <f>' and as shown in the 
Eq. ( p.9|) . Thus we can schematically indicate the azimuthal dependence of the potential 
matrix elements as 

V]ff (q ; , q) = VZf {e lJU( ^, cos(0' - 0)} (3.12) 
For the special case q = z this reduces to the simpler form 

VfM, q) = e^'^V^(q', q, 9'), (3.13) 
which is the driving term in the coupled set of Eqs. (|2.21|) . 



S 



We assume that this simple dependence on cf>' and <fi given in ( |3.13j ) carries over to the 
solution of the integral equation and choose as ansatz 

Tfl (q' ; q) = e*W-*)TZ*tf, q, 9') (3.14) 

Inserting this into Eq. Q2.15| ) one obtains 

T^, q z) = e^'-^iy, q ,e') 

+7 £ / dY Vg£ {e iM ^"\cos{<t>' - 0"), </, q"} 

4 A" J 

G (q'y^''-^i(q'',q,9'') 

j poo rir p2iT 

+ V AW, </,) E/ dq" q" 2 d9" sin 9" d(f>" 
4 A „ Jo Jo Jo 

Vfl {e iM ^"\ cos(0' - 0"), g', g"} G (g") 

e iA W'-^T${{q",q,e") (3.15) 

The integrand is periodical with respect to <fi", with the period being 2n. Consequently one 
can set 0' = for the ^"-integration. This leads to 



T^q', qz) = e 



V^(q',q,9') 



" ' jn 1 d9" sin 9" 
o Jo 



2 71 



Vf A * {e-^",cos 0", q', q"} Go{q")e^" T^l{q\ q, 9") 



(3.16) 



verifying the correctness of the ansatz of Eq. (|3.14j ). 
Inserting this result into Eq. Q2.21|) gives 

Tft(q',q J 9') = VZ s A t (q',q,9>) 

+ o / <W' 2 / d(cos6") 
I Jo J-l 

vl?t' A (q', q", 9', 9")G (q")T^(q", q, 9") 

I roc rl 

- / dq"q" 2 / d(cos0") 
4 jo J-i 

vl%*tf, q", 9', 9")G (q")T T(q", q, 9") 



+ 



with 



/•27T 

vZf&tf, q", 9', 9") = / d^'e-^'-^V^M, q") 



(3.17) 



(3.18) 



The driving term of Eq. (|3.17|) is a special case of Eq. ( |3.18|) for 9" = and A" = A (up 
to a factor (2n)~ 1 ). In summary, for the case S — 1 we end up with two coupled integral 
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equations in two variables for given A- values (only A = 1,0 are necessary). The azimuthal 
integration over (j)" can be performed independently, and does not enter the integral kernel. 

The case S = is much simpler, one has only one integral equation in two variables, 
similar in its structure to the two-boson case discussed in Ref. |§. 



IV. REPRESENTATION OF NN POTENTIALS 



In this Section we want to demonstrate the relative ease with which modern NN poten- 
tials that are given in operator form can be incorporated in a three-dimensional formalism. 
Our choices of NN potentials are a Bonn one-boson-exchange (OBE) potential Jl9| in the 



parameterization of Bonn-B JT5| and the Argonne coordinate space potential AV18 ||16|| . 

The OBE potential consists of pseudoscalar-, scalar-, and vector meson exchanges, de- 
rived from the corresponding Feynman diagrams. A three-dimensional reduction of the 
Bethe-Salpeter equation is as achieved via the Blankenbecler-Sugar reduction. Details can 



be found in Refs. [fT5,|19|. The resulting potential operators are given as 



nq',q) 
^(q',q) 

nq',q) 



.</: 



9 



(2tt 

2 



u(q')u(q)u(— q')w(— q x ^ ' 



q) 



ps 



(2tt) 3 

^[(q'-q) 2 



q' - q) 2 + A 

— u(q )Yu{q)u(-q )7°u(-q) 7-7 
E q' 



q) 2 + /i 2 



ps 



q' 



q) 2 + H 2 V 



E (2tt) 



S-Xq')7 M w(q)w(-q')7M u (-q) 



(4.1) 
(4.2) 
(4.3) 



+ 



{4m 2 M(q')7 M M(q)M(-q')7MM(-q) 

-2m«(q , )7^(q)«(-q')[(^ - £)G# ~ 
-2mu(q')[(E' - E){g^ - 7 M 7°) + (ft 



7 M 7 



+ (P2+P 2 )^]u(-(l) 

p 1 r}u(q)u(-q')-f tl u(-q) 
+u(q')[(E' - E){g^ - 7^7°) + fa +p£)>(q) 
«(-q')[(£' - £)(<?J - 7,7°) + (P2 +p' 2 )>(-q)} 

9vfv 



2m 



7m7 ) + (P2+^) w 

{4m n(q')7 M ^(q)^(-q')7M w ( _ q) 



-qj 



-«(q')7^(q)«(-q')[(^ - ^ - 7.7°) + (ft + p' 2 )/J 
- «(q')[(^ - W - 7 M 7°) + (ft +Pi)1«(q)«(-q')7M«(-q)} 

Here m stands for the nucleon mass. In the case of the vector potential one has (ft +p l 1 ) tl = 
(E + E', q + q') and (ft + P2T = (E + E', — q — q'). The coupling constants g ps , s ,v and f v , 
the cutoff functions F and the meson masses are given in fli"5| . 

In order to bring this OBE potential in a form consistent with our three-dimensional 
(3D) equations, the bilinear Dirac forms have to be expressed in terms of the operators 
Wi of Eqs. QjUp . The result in form of the Wolfenstein operators of Eq. ( |3.5p is given in 
Appendix B. 

Another often used modern NN potential is the Argonne potential AV18 ||16|1 . It is 
originally presented in configuration space and has the general form 
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V(NN) = v EM {NN) + v*(NN) + v R (NN). 



(4.4) 



Here v EM (NN) represents an electromagnetic part, which we omit in this work. The one- 
pion-exchange (OPE) part v^lNN) is charge dependent and has the standard form. The 
Yukawa and tensor functions contain exponential cutoffs, thus do not have an analytical 
Fourier transform to momentum space. The intermediate- and short-range phenomenological 
part v R (NN) is expressed as a sum of central, tensor, spin-orbit, L 2 , and quadratic spin- 
orbit pieces (abbreviated as c, t, Is, 12, ls2, respectively) in different spin (S) and iso-spin (T) 
states: 

v§ T (NN) = v c ST (r)l + 4t(05i 2 + v l i T (r)L • S + v l2 T (r)L 2 + v l s s2 (r)(L • S) 2 , (4.5) 

where Si 2 denotes the standard tensor operator. The specific form of the radial functions 
as well as the potential parameters are given in Ref. |p.G(| . For applying this potential in 
our formulation, we need to perform the transition to momentum space. For the terms 
contributing to v R (NN) we obtain explicitly 



4 s t (q',q) 



vir(q',q) 



i 



2tt 2 Jo 
1 ( 



drr 2 j (pr) v c ST (r) 
-3 (T 1 ■ (q' - q) <r 2 • (q' - q) 



2tt 2 

% 

2^ 
1 



+ a- 1 ■ <x 2 



drr 2 j 2 (pr)v t ST (r) 



1 f°° 

S-(qxq')- drr 3 ji(pr)t4 s T ( r ) 
p Jo 

,2 



Tr^q ■ q , 

2n z p Jo 



2tt 2 



12 2/1 2\ 

q q (i - 7 ) 



drr ji(pr)v ST (r 
1 



p 2 Jo 



drr j2(pr)v ST (r) 



1 / \ 1 z* 00 

^l(q', q) = ■ q " S ■ qS • q') - ^ drr 8 *^)^) 



2tt 2 V2 



1 



Triq x q 



/\2 



-o-i • (q X q)cr 2 ■ (q x q 



1 



p 2 Jo 



/ drr A j 2 (pr)v l s S T(r) 



(4.6) 
(4.7) 
(4.8) 



(4.9) 



(4.10) 



where p = |q' — q|. The resulting operators can be easily represented as function of the 
Wolfenstein operators W%. The final expression are given in Appendix C. 



V. CONNECTION TO A PARTIAL WAVE REPRESENTATION 

In order to compare with standard partial representations, especially NN phase shifts, 
we need to make connection to the standard partial wave representation. The partial wave 
projected T-matrix element is defined as 

Tff\q) = (q(l'S)jmtm t \T \q(lS)jmtm t ) , (5.1) 
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where the states \q{lS)jmtm t ) are given as 

\q{lS)jmtm t ) = E C(lSj, \im — jj) \qlfi) \Sm — n) \tm t ) . (5.2) 

We choose the standard normalization for those states, namely 

(q'(l'S')fm't'm' t \q(lS)jmtm t ) = g ^ Siq6 s , s 5 jrj 6 m , m 5 t , t 5 m > tmt . (5.3) 

The connection to the states |q; qSA) can be found using Eq. (|2.9|), and one finds 
\q;qSA) = J2 WS)jm) 

Ijm 

E C(lSj; /im - ^(q)e^ m -^d s m ^ A (9). (5.4) 

Consequently one obtains for the T-matrix 

7X5(q', q) = V2(l - ^(-) S+t ) (q'| q'SA'\ (*| T \t) |q> |$SA> 

7T 

= ^(1 " VnH S+t ) E + %(-)")(! + %(-)') 

Z'Zj'm 

E C7(J'£j; //, m - //) ^(^^'^-^a^') 
i*' 

E C(^j; /i, m - /u)Y^ e -^-^4_ M (0). (5.5) 

For q parallel to the z-axis, i.e. q — z, one finds after some straightforward algebra the 
on-shell relation 

Tg£{q,q,<n = ki-vA-) s+t )Y, T iTm+vA-) l ')(i+v«n l ) 

/ I'lj 

fop _|_ 1 /?/ 4- 1 

- -C7(Z'Sj; OA')^,^ ^—C(lSj; OA). (5.6) 



Here the relation 



*V(<f) = f^KW^) = ^^'(-r'D^e'O) (5.7) 

was used together with an addition theorem for D-functions. 

Once the partial wave projected T-matrices are calculated, we can use the the well known 
connection to the partial wave S-matrices 

S$( t = 8 vl -™rn4I$i\ (5.8) 

which are parameterized by the standard partial wave phase shifts. 

The inversion of Eq. (|5.6|) can be accomplished by using the orthogonality relation 
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^ # jT ^ sin ^(^0)2^(^0) = ^7<W ( 5 - 9 ) 



In addition we take into account Eq. ( 3.14 ) for the on-shell T-matrix and find after a few 
steps 

T?,i\q) = 4tt(1 - ^(-) s+t )- 1 (l + + Hr(-)')- 1 

x/ ^?/f ITI E 0A')C(Z^; OA) 

4? + 1 A' A 

dicose^i^Tf^q^'). (5.10) 



Finally we connect the physical T-matrix element to the partial wave projected S- 
matrices. We insert (5.6) into ( 2.2SQ . In proceeding we use 

Di* A ,(4>'e'0) = n A - A 'e lA *'d^A0>) (5.11) 



and an addition theorem of D-functions together with (|5.7f) . Then we end up with the 
standard form 

a ( v x v 2 m[ m' 2 qq' \ T \ v x v 2 m x m 2 q) a 

2 y47T Ttimq ^22 22 22 



V27+lC(lSj, 0Ao)C(l'Sj, A - A' 0! K)Yi>,a -a' </>') (5.12) 

One can choose 0' = by putting q' into the x — z plane 

In order to unambiguously define the normalization we also give the expression for the 
spin-averaged differential cross section for nucleon species v%v 2 as 

^ = (2tt) 4 [^j i E \a{^2rn' J ^m' 2 q^\T\v 1 v 2 m 1 rn 2 (^ a \ 2 (5.13) 



VI. RESULTS AND DISCUSSION 

In order to demonstrate the feasibility of our formulation when applied to NN scattering, 
we present numerical results for two NN potentials of quite different character, the Bonn-B 



15 1 and the AV18 |16] potential models. For calculating the NN t-matrix we solve the 



integral equation, Eq. ( 3.171) . For S=0 this is a single, two-dimensional integral equation in 



two variables, for the parity even and odd part, respectively. For S=l there are two sets of 
coupled integral equations in two variables for the helicities A = 0,1. The 0"- integration 
over the potential, Eq. (|3.18|), can be carried out independently, and thus does not enter the 



integral equation as separate variable. As it turns out, both potentials only depend weakly 
on this angle, and one typically needs 10 integration points in the case of Bonn-B and 16 
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in the case of AV18 to have sufficient accuracy. Typical grids needed for the ^"-integration 
consist of 32 to 48 points and for the (/'-integration of 48 to 72 points, depending on the 
desired accuracy. 

The Cauchy singularity in Eq. ( |3.17 ) is separated into a principal value part and a 5- 



function part, and the principal value singularity is treated by subtraction. In case of the 
Bonn-B potential the integration interval for the g"-integration is covered by mapping the 
Gauss- Legendre points u from the interval (0, 1) via 

q = 6tan f^^J (6-1) 

to the interval (0, oo). Typical values of b are 1000 MeV/c. For the AV18 potential this 
type of integration map is less useful due to the difficulty in handling the numerical Fourier 
transform for very large values of q", and the integration interval is cut off at 150 fm _1 , 
which is according to our experience sufficient. The Fourier-Bessel transform of the AVI 8 



potential functions are carried out using Filon's quadrature formula [21], which proves to 
be more accurate in handling the strong oscillatory behavior of the integrands in Eqs. ( f4.6p - 
( |4.10| ) for large values of the integration parameter r, compared to e.g. Simpson's rule. 

A stringent test for our numerics is a comparison of the NN phase shifts, which we 
obtain from Eq. (|5.10|) together with Eq. ( |5.8j ), with phase shifts calculated with standard 
partial wave techniques. The result for the Bonn-B phase shifts is given in Table I for 
projectile energies 100 and 300 MeV. The agreement of both calculations is excellent. In 
Table II the equivalent comparison is given between the phase shifts calculated from the 
AV18 potential with our 3D formulation and a standard partial wave calculation (in this 
case in coordinate space). The agreement is also very good, however not as perfect as in 
the Bonn-B case. The reason for this slight discrepancy is presumably twofold. First, we 
carry out a numerical Fourier-Bessel transformation of the different potential functions of 
AV18. The solution of the t-matrix integral equation requires an evaluation of the potential 
function on a grid of the size n^n x ngn x n q n x n q i x ngi. For computational economy 
we calculate the potential functions on a fine grid for p = |q" — q| and obtain the points 
actually needed in the calculation via interpolation. This procedure naturally leads to larger 
numerical errors compared to a direct evaluation of algebraic expressions as is the case for 
Bonn-B. The differences can be clearly seen when comparing Tables I and II. In both cases 
a comparable grid for the t-matrices is used. 

Next we want to display the angular and momentum dependence of the half-shell t- 
matrices TJ[^(q, q , 6) as generated from Eq. ( |3.17| ). For the case S = this is done in Fig. 1 
for the Bonn-B potential for the parity-even (left side) and for the parity-odd case (right 
side). The parity-even t-matrix exhibits a similar behavior to the symmetrized t-matrix of 
the two-boson case studied in Ref. ||. It shows strong forward and backward peaks around 
the on-shell momentum qo. Once far away from go the momentum dependence is quite 
moderate. The parity-odd t-matrix also displays the the strong peaks around the on-shell 
momentum, however the forward and backward peaks have opposite signs. 



For the case S = 1 there are two coupled integral equations as shown in Eq. ( |3.17 ). Due 



to the symmetries of the potential, Eq. ( |2.19p , it is sufficient to consider only the helicity 
combinations A', A = 0,1. The corresponding half-shell t-matrices, T^(q,qo,9) obtained 
from the Bonn-B potential are displayed in Figs. 2 through 5. The real and imaginary 
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parts show a rich structure as function of angle and for momenta smaller than 1000 MeV/c. 
For momenta larger than 1000 MeV/c all amplitudes show no or only a very weak angular 
dependence. We would like to point out that at 300 MeV projectile energy partial waves have 
to be summed at least to J = 6 to obtain a converged result for calculating e.g. observables. 

A further question of interest is how strong the difference is between the half-shell t- 
matrix amplitudes derived from the two potentials presented here. In Fig. 6 the real and 
imaginary parts of TJ[, s ^(q, qo, 6) for S — calculated with AV18 are shown. Of course along 
the line q = qo the corresponding amplitudes derived from the two potentials are identical. 
However, if one compares the detailed structures of especially the parity odd amplitudes, 
one sees that at large momenta q the Bonn-B amplitude shows more structure than the 
one derived from AV18. For the case S=l we only want to show a selection of t-matrix 
amplitudes calculated from AV18. In general one can say that some amplitudes are very 
similar to the ones derived from Bonn-B, others are quite different. The choices displayed 
in Fig. 7 are being representative. The upper two t-matrices show a considerable difference 
to the corresponding ones of Bonn-B, while the lower ones are relatively similar. 

Next we turn to the on-the-energy-shell t-matrix amplitudes as given in 
Eq. ( |2.28|) in terms of the T^f^ quantities. Here it is interesting to consider 
\ a (^i^2 m i^2QoQ\T\i / i^2 ri T'i ri T'2 < io)a\ 2 = I (m^m^ \T\rri1m2) | 2 , where m = ±|. If one takes ro- 
tational symmetry and parity invariance into account, one ends up with six independent 
amplitudes, 
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(6.2) 

The six amplitudes given above are displayed in Figs. 8 to 10 for Bonn-B as function of 
energy and cm. scattering angle cos#. Though both potentials are only meant to be 
applied below the pion production threshold, we show the on-shell amplitudes up to 1 GeV. 
Since we do not work with partial waves, a calculation at higher energies takes the same 
effort as one at very low energies. The sum over all indices of the on-shell amplitudes 
shown gives the spin averaged differential cross section as indicated in Eq. Q5.13p . We would 
like to remark that these on-shell amplitudes can also be obtained from the partial wave 
projected S-matrix elements as indicated in Eq. ( |5.12j ). We used this relation for numerical 
tests of our formulation. 



VII. SUMMARY 

Two nucleon scattering at intermediate energies of a few hundred MeV requires quite a 
few angular momentum states in order to achieve convergence of e.g. scattering observables. 
This is even more true for the scattering of three or more nucleons upon each other. An 
alternative approach to the conventional one, which is based on angular momentum decom- 
position, is to work directly with momentum vectors, specifically with the magnitudes of 
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the momenta and the angles between them. We formulated and numerically illustrated this 
alternative approach for the case of NN scattering using two realistic interaction models. 
The momentum vectors enter directly into the scattering equation, and the total spin of the 
two nucleons is treated in a helicity representation with respect to the relative momenta of 
the two nucleons. Using rotational and parity invariance one finds in the triplet case S=l 
a set of two coupled Lippmann-Schwinger equations for each parity and each initial helicity 
state. Because of symmetry properties only two of the originally three initial helicity states 
need to be considered. In the singlet case S=0 there is only one single Lippmann-Schwinger 
equation for each parity. The Lippmann-Schwinger equations (uncoupled and coupled) are 
two-dimensional integral equations in two variables for the half-shell t-matrix and three 
variables for the fully-off-shell t-matrix, namely two magnitudes of momenta and one an- 
gle. Though we start with a three-dimensional integration, one angle can be integrated out 
independently, so that we are left with a two-dimensional integral equation. 

A formulation without angular momentum decomposition is best suited for interaction 
models which are given in an operator form. In this work we considered the Bonn-B and the 
AV18 potentials. The helicity representation of the potentials is given in the six linear 
operators Qi of Eq. ( |3.6| ), depending on spin and momenta, which are most suited for 
our formulation. For completeness we also give the potentials represented via Wolfenstein 
operators. We demonstrate the feasibility and accuracy of our three-dimensional formulation 
by projecting the on-shell t-matrix elements on angular momentum states and compare the 
resulting phase-shifts with those obtained from standard partial wave projected Lippmann- 
Schwinger equations. The agreement is very good and demonstrates the numerical reliability 
and accuracy of our method. 

The complete set of half-shell t-matrix amplitudes at 300 MeV laboratory energy is 
displayed for the Bonn-B potential to give an impression of their rich structure in momentum 
and angle variables. Nevertheless the resulting angular dependence is smoother than the 
individual contributions resulting from the higher angular momentum states. 

In some representative examples it is shown that those half-shell amplitudes can differ 
significantly depending on the potential employed. Of course, the observable quantities, the 
on- shell amplitudes, are equal since they describe the same NN data base. We display the 
magnitudes of the physical t-matrix elements derived from the Bonn-B model as a function 
of scattering angle and energy. Though the NN potential is tuned only up to about 350 
MeV our results are shown up to 1 GeV to demonstrate the ease by which also high energies 
can be handled in our approach without partial waves. 

We would like to emphasize that the here developed scheme is algebraically very simple 
to handle provided potentials are given in an operator form. This is e.g. the case for all 
interactions developed within a field theoretic frame work. In addition, the solution of 
two-dimensional integral equations does not pose any difficulty for modern computers. 

This work is intended to serve as starting point towards treating three-nucleon scattering 
without partial waves. In Refs. ||[7j three-nucleon bound and scattering states have already 
been treated without partial waves, however neglecting spin degrees of freedom. In the 
present work we show at the two-nucleon level that the inclusion of spin degrees of freedom 
is quite readily possible. Future work needs to combine the experience with the NN system 
gained here with the three-body calculations without partial waves of Refs. ||[7j . This will be 
an important step forward since present 3N scattering codes based on partial wave methods 
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reach the limit of todays supercomputers already at intermediate energies below the pion 
threshold. Future applications at higher energies require new techniques. 
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APPENDIX A: THE TRANSFORMATION MATRIX A 



In this appendix we explicitly give the transformation from the Wolfenstein operators 
given in Eqs. (|3.1| ) to the operators of Eqs. ( |3.6| ), i.e. the matrix indicated in Eq. (|3.7|). 
The matrix is explicitly given as 



A 
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(Al) 



The matrix elements Aij are scalar functions. We calculate 



detA 



qq 



7\/i - i 2 Vlq' + q||q' 



(A2) 



and find the inverse matrix, A-,- to be 
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(A3) 



For both matrices we introduced the abbreviations 



a 
b 

c 
e 

f 



1 - 7 2 
q' 2 + q 2 - 2q'q~f 
q' 2 + q 2 + 2^7 



(g' 2 + g 2 ) 2 - 4g /2 g 2 7 2 
I 1 
g' 2 q 2 



9=^2--^ 



APPENDIX B: THE OBE POTENTIAL AS FUNCTION OF WOLFENSTEIN 

OPERATORS 

In terms of the Wolfenstein operators Wi of Eq. (|3.1|) , the Bonn-B potential takes the 
following form: 
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v ( ' \ g2ps [™ !™ F ^ - q ^ ° ps 

Mq ' qj (27r) 3 4m 2 V E'V E (q' - q) 2 + m 2 ps WW { ' 

vM n] _A /™ /^ *?[(q'-q) 2 ] a m9 x 

= i rw rm F 2 [( q ' - g) 2 ] (g 2 Av + 2&,/ t ,g t + / 2 o tf ) 

" iq ' qj (27r) 3 4m 2 V£'V£(q'-q) 2 + m 2 WW ' 1 ' 

where m denotes the nucleoli mass and W = \/m 2 + q 2 . The operators ps>s ^ v are given as 

PS = -\[{W - W) 2 {q' 2 + q 2 + 2q' qi )W 4 + (W + W) 2 (q' 2 + q 2 - 2g'g 7 )W 5 ] 

~(W' 2 - W 2 )^(q' 2 + g 2 ) 2 - 4g' 2 g 2 7 2 W 6 (B4) 

S = -(WW - g'g 7 ) 2 Wx - i{W'W - q'qi)q'q\Jl - 7 2 W 2 

+g'V(l-7 2 )W 3 (B5) 
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For the actual calculations the above given terms are re-expressed in terms of the operators 
Qi via Qj = X)j A^Wj an d inserted into the L.S. equation Eq. (|3.17|) . These reformulations 
are straightforward, and can e.g. be carried out via symbolic manipulation packages. 

APPENDIX C: THE ARGONNE AV18 POTENTIAL AS FUNCTION OF 

WOLFENSTEIN OPERATORS 



In terms of the Wolfenstein operators W,i of Eq. ( |3.1| ), the AV18 potential takes the 
following form: 
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In our actual calculation the potential enters described by the set of operators f2j of Eq. Q3.6|) . 
In order to display the relative simple form the potential takes in these operators, we also 
want to explicitly show them. 
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TABLES 



TABLE I. Comparison of the nucleon-nucleon phase shifts calculated from our 
three-dimensional formulation (helicitiy) with a standard partial wave calculation (carried out 
in momentum space) for the Bonn-B potential at 100 and 300 MeV laboratory energy. 





E lab = 100 MeV 


Em, = 300 MeV 


2S+1 L 


3D helicity 


partial wave 


3D helicity. 


partial wave 




25.1928 


25.1929 


-8.1755 


-8.1756 


1 u 


9.8046 


9.8046 


-11.4799 


-11.4799 


1 L 


-16.3131 


-16.3451 


-28.6946 


-28.8747 


3 Pi 

± L 


-13.4677 


-13.4677 


-26.3800 


-26.3800 




41.9858 


41.9870 


4.0667 


4.0676 


-LSI 


-12.9847 


-12.9846 


-23.7182 


-23.7181 


° 1 


-2.2360 


-2.2357 


-4.0268 


-4.0265 


1 J D2 


3.3411 


3.3411 


7.4888 


7.4888 


3 Do 


17.6710 


17.6710 


25.3616 


25.3617 


3 Po 


11.7356 


11.7356 


17.3981 


17.3981 


3 Fo 

1 2 


0.7705 


0.7705 


0.5236 


0.5238 


Fn 


2.8402 


2.8402 


2.0166 


2.0166 


^3 


-2.4397 


-2.4397 


-5.5865 


-5.5865 


3 F 3 


-1.6484 


-1.6484 


-4.0097 


-4.0097 


3 D 3 


0.4203 


0.4855 


2.5719 


2.5720 


3 G 3 


-1.0105 


-1.0105 


-4.4051 


-4.4051 


£3 


-3.6604 


-3.6604 


-7.2233 


-7.2233 


'G 4 


0.4092 


0.4092 


1.3556 


1.3556 


3 G 4 


2.2624 


2.2624 


7.3000 


7.3000 


3 F 4 


0.4203 


0.4203 


2.4491 


2.4491 


3 H 4 


0.1082 


0.1082 


0.5077 


0.5077 


e 4 


0.5575 


0.5575 


1.5509 


1.5509 
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TABLE II. Comparison of the nucleon-nucleon phase shifts calculated from our 
three-dimensional formulation (helicity) with a standard partial wave calculation (carried out in 
coordinate space \22\) for the AV18 potential at 100 and 300 MeV laboratory energy. 





E lab = W0MeV 


E lab = 300 MeV 




3D helicity 


partial wave 


3D helicity. 


partial wave 


So 


25.99 


25.94 


-4.62 


-4.60 


3 Po 


8.69 


8.69 


-11.05 


-11.06 


'Pi 


-14.19 


-14.20 


-26.18 


-26.28 


3 P1 


-13.06 


-13.07 


-28.38 


-28.49 


3 s, 


43.69 


43.56 


8.15 


8.16 




-12.08 


-12.09 


-24.80 


-24.90 


£1 


-2.49 


-2.49 


-4.38 


-4.39 


X D 2 


3.81 


3.81 


9.45 


9.44 


3 D 2 


17.14 


17.10 


25.11 


25.02 


3 P 2 


11.02 


11.00 


16.96 


16.91 


3 F2 


0.67 


0.67 


0.77 


0.76 


£2 


2.70 


2.70 


2.21 


2.21 


^3 


-2.23 


-2.23 


-4.87 


-4.88 


3 Fs 


-1.35 


-1.35 


-2.51 


-2.51 


3 D 3 


1.61 


1.61 


5.22 


5.21 


3 G 3 


-0.93 


-0.93 


-4.19 


-4.20 


£3 


-3.50 


-3.50 


-7.17 


-7.16 




0.40 


0.40 


1.42 


1.42 




2.22 


2.22 


7.35 


7.34 


3 Fa 


0.45 


0.45 


2.75 


2.74 


3 H 4 


0.07 


0.07 


0.31 


0.31 


e 4 


0.51 


0.51 


1.54 


1.54 
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FIGURES 



FIG. 1. The half-shell t-matrix T][%(q,qo,6) for q =375 MeV/c, corresponding to 
Ei a b=300 MeV, calculated from the Bonn-B potential. The left side displays the real and imaginary 
parts of T 1 q 1 (q, <7oj &) as function of q and x = cos 9. The right side shows T^ 100 (q, q , 6) . The 
units of Tffi (q,q ,0) are 10" 7 MeV" 2 . 

FIG. 2. The real part of the parity-even half-shell t-matrix T^}^(q,qo,6) for qo=375 MeV/c. 
The helicities A', A = 0,1 are kept as subscripts of T. All t-matrices are calculated using the 
Bonn-B potential and are given in units 10 -7 MeV -2 . 

FIG. 3. Same as Fig. 2, but for the imaginary part of T^(q,qo,6). 

FIG. 4. The real part of the half-shell parity-odd half-shell t-matrix T^ n (q, qo, 6) for 
go=375 MeV/c calculated from the Bonn-B potential. The description is the same as in Fig. 2. 

FIG. 5. Same as Fig. 4, but for the imaginary part of T'^}^ 1 {q,qQ,6). 

FIG. 6. The half-shell t-matrix T$%(q,q ,6) for q =375 MeV/c calculated from the AV18 
potential. The left side displays the real and imaginary parts of T^ 1 (q,qo,9) as function of q and 
x = cos6. The right side shows T^ 100 (q, q , 6). The units of T^ t {q l , qo, 8) are 1(T 7 MeV -2 . 

FIG. 7. Selected half-shell t-matrices T£}£(q,qo,0) for q =375 MeV/c calculated from the 
AV18 potential. The helicities A', A = 0, 1 are kept as subscripts of T, which is given in units 
of 1(T 7 MeV -2 . 

FIG. 8. The absolute value squared of the physical on-shell T-matrix elements 
\{m' 1 m' 2 \T\mim2)\ 2 in units 1CT 14 MeV -4 . The values of m = ±\ are abbreviated as ±. 

FIG. 9. Same as Fig. 8, but for different m combinations. 

FIG. 10. Same as Fig. 8, but for different m combinations. 
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